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Abstract 

We investigate the 2-point correlation function for the vector current. The 
gluons provide dressings for both the quark self energy as well as the vec- 
tor vertex function, which are described consistently by the rainbow Dyson- 
Schwinger equation and the inhomogeneous ladder Bethe-Salpeter equation. 
The form of the gluon propagator at low momenta is modeled by a 2-parameter 
ansatz fitting the weak pion decay constant /tt. The quarks are confined in the 
sense that the quark propagator does not have a pole at timelike momenta. 
We determine the ground state mass ttlq in the vector channel from the Eu- 
clidean time Fourier transform of the correlator, which falls off as e"™"^ at 
large T. mo lies around 590MeV and is almost independent of the model form 
for the gluon propagator. This method allows us to stay in Euclidean space 
and to avoid analytic continuation of the quark or gluon propagators into the 
timelike region. 

ll.lO.St, 12.38.Aw, 14.40.Cs 
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I. INTRODUCTION 



The purpose of the present work is the determination of the ground state mass of the 
vector qq bound state within the Dyson-Schwinger (DS) and Bethe-Salpeter (BS) approach 
|Jl|-P]. The starting point is the 2-point correlation function 

^^.uiq) = I f/W«-"(TJ^(x)J,(0)) (1) 

for the vector quark currents J^{x) = q{x)'y^q{x), with q{x) a quark field. Within the ap- 
proximations used in the present paper, which will be described below, U^i, is determined by 
the dressed quark propagator {Tq{x)q{0)) and the vertex 3 point function {Tq{x)J^{y)q{0)). 
The objective of this work is to study the large-time behavior of the 2-point correlator in 
Euclidean space in order to determine the mass of the lowest vector meson, as in lattice 
gauge calculations ^ and instanton models ^ . QCD sum rules also consider the Euchdean 
2-point correlator, but attempt to determine the ground state mass by using the Borel trans- 
form to isolate the lowest mass pole |^ . Our method should be seen in contrast to an earlier 
approach, which solves the BS equation on shell 0. 

We now briefly review the DS-BS approach used here. The nonperturbative nature of 
QCD at low and intermediate momentum transfers makes the use of effective models a 
natural and necessary tool. During the past few years strong progress has been made in 
the framework of the Dyson-Schwinger (DS) approach which is based on a coupled 

set of relations between dressed quark, gluon and ghost propagators and vertex functions. 
In order to handle this system it is of course necessary to make certain simplifications and 
truncations. One approach, which is commonly referred to as rainbow approximation, uses 
an undressed quark-gluon vertex and solves the Dyson-Schwinger (DSE) for the dressed 
quark propagator G{p) with a given dressed gluon propagator D^^^^{p) as input (c.f. Fig.[^): 

= / -0^,9s'D;lip - g)7.f G(g)7.| • (2) 
The quark propagator G and the quark self energy E are related by: 

G(p)-' = t^.p + j:(p) . (3) 

The form of the gluon propagator D^^^ij)) for small momenta is basically unknown from 
QCD and therefore has to be modeled while following certain requirements which will be 
discussed below. In case of this special truncation it is possible to formulate an effective 
field theory using functional integral techniques, which is known as the Global Color Model 
(GCM) This has the striking advantage that in the chiral limit (zero current mass 

for the light quarks, i.e. Mu = Md = Mg = 0) one is able to establish an easy connection 
between the dressed quark self energy leading to dynamical chiral symmetry breaking on 
one side and the existence of massless Goldstone bosons as pseudoscalar qq bound states, 
on the other side. As a consequence it is possible to perform a systematic chiral low energy 
expansion in terms of composite Goldstone boson field and therefore make a connection 
between the phenomenological chiral hadronic theory (Chiral Perturbation Theory) |jlO|] and 
an effective quark based theory of QCD, the GCM. 
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Furthermore one is also able to study the coupling of composite qq systems to an external 
electromagnetic gauge field while guaranteeing local U{1) gauge invariance which is manifest 
in the corresponding Ward-Takahashi identity [[I^jH- It turns out that the chiral invariance 
and electromagnetic current conservation, which are both essential for studying low energy 
hadronic phenomena, are very difficult to maintain if one tries to go beyond the rainbow 
approximation and are lost in general if model ansatze for the quark-gluon vertex are 
used. 

Though the rainbow approximation and the GCM naturally violate local color SU{3)c 
gauge invariance and renormalizability, they provide a very successful description of vari- 
ous nonperturbative aspects of strong interaction physics and the QCD vacuum as well as 
hadronic phenomena at low energies. These include for instance quark confinement 0, QCD 
vacuum condensates [p|,p!3|-pr5[] , Ua{^) breaking and the rj — r]' splitting |T6[, low energy chiral 



dynamics of Goldstone bosons (tt, K, rj) |[T7|j9|,pr^- pO| , p!5| , meson form factors heavy-light 



mesons systems at finite temperature or soliton |^ and Fadeev descriptions 
of the nucleon. 

The method is more involved when applied to higher mass states instead of Goldstone 
bosons, e.g. p, u, a or A mesons. A mesonic qq bound state with mass m in a channel with 
spin-fiavor index 6 is determined in the ladder approximation as solution of the homogeneous 

fi,(P, q' = -m') = J ^,9s'D;tiP - K)^,^G{K^)^e{K, q' = -m')GiK_h,^ , 

(4) 

where we have used the notation K± = K The ^e{P, q) is related to the Bethe-Salpeter 
amplitude Xe{P,<l) as: 

Xe{P, q' = -m') = G{P+)Qe{P, q' = -m')G{P^) . (5) 

As mentioned above in the pseudoscalar channel the existence of a massless state or in 
other words a solution of (^) at = is automatically guaranteed in the chiral limit if 
the rainbow DSE is satisfied. For a small finite current quark mass Mq a perturbative 
expansion in Mq and can be performed. This is however not possible for larger meson 
masses and therefore an explicit numerical solution of the corresponding homogeneous BSE 
(^) is required. The trouble hereby is that the meson momentum is defined in the timelike 
region, = — m^, which in turn requires the dressed quark propagator G{p) in the timelike 
region. However G{p) is a priori only determined at spacelike p from the DSE (|^), because 
the model gluon propagator D{p) is defined only for spacelike Euclidean momenta. It is 
therefore necessary to perform an analytic continuation of the dressed quark propagator 
G{p) from Euchdean space into Minkowski space p4 — * ipo [0]. This is a very dangerous and 
unsafe procedure because it requires in fact the knowledge of the singularity structure of 
G{p) in the whole complex plane, which we do not have. In Ref. an alternative method 
has been proposed which stays completely in Euclidean space avoiding analytic continuations 
of any of the dressed propagators. It is the analog to what is used in instanton models 
and lattice QCD calculations |Q where working in Euclidean space is essential. 

The 2-point correlation function 
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^e,eM = J d^xe''^-^{TJe,{x)JeM) (6) 

with quark currents Je = qTeq, where 9 stands for the spin-flavor of the operator Tq, is 
determined as a quark loop containing the dressed vertex function rg{P,q) (c.f. Fig.^ for 
Te = The Tg^P, q) itself is given as the solution of the inhomogeneous ladder BSE (c.f. 
Fig-I): 

r,(P, q) = -zTe + J ^,9s'D^^AP - Kh,^G{K^)Te{K, q)G{K.)^^^ (7) 

with the inhomogeneity —iTg. This form can also be written as 

Te{P,q) = -iTe + Te{P,qf'', (8) 

where Tq is the bare vertex and rg{P,q)^^ is the nonperturbative dressed vertex. The 
nonperturbative parts of the vertex functions can be constrained by the forms of the nonlocal 
condensates which also has been used in our recent work Moreover, it has been 
shown that this formalism can be used for the nonperturbative part of the vector vertex 
function used at low momentum transfer for nuclear magnetic dipole moments [^. This 
gives possible constraints for the dressed vertex function treated in the present work, which 
we discuss below. 

The momentum q is in general off shell and can be either space or timelike. Close 
to the mass shell q^ ~ — the rg{P,q) and the solution of the homogeneous BSE 
fle{P, q^ = —m'^) are related by 

Ve{P,q^^-m^)^ 7^ , '- . 9 

qz _|_ 

In our study the q^ will be always spacelike. Other than the dressed quark or gluon propa- 
gator the correlator 'n.Q^0^{q^) has a known Kallen-Lehmann spectral representation Its 
poles are at q^ = —mi'^, I = 0, 1, . . ., where the mj are the meson resonance masses and 
the pole residua are basically given by the couplings fi between the interpolating current Je 
and the on-shell states |J). Otherwise Uoj^e^il'^) is holomorph in the complex plane and falls 
off sufficiently fast for oo, so that Cauchy's theorem can be applied. A convenient 

method to filter out the lowest state with mass mg is to take the the one dimensional Fourier 
transform of with respect to the time component of the Euclidean momentum q^ 

n;,,,(T) = / ^e^^^^Ue.oM, Q = 0) (10) 

which falls off as e"'""^ for large Euclidean times T. 

It is obvious that this method will work in reality only if the ground state mass mg 
is sufficiently well separated from the mass of the first excited state mi. If mo and mi 
are too close it will be numerically difficult to extract the exponential falloff, because this 
situation will require numerical very large values for T in order to see the falloff and for those 
large values the signal is likely to be swamped by numerical errors due to highly oscillating 
behavior of e*''""^. On the other hand it is clear that we have indeed avoided to make any 
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assumption about the analyticity of the propagators. Instead we are using the analyticity 
structure of the 2-point correlator Uo-^o^lq'^) in the complex plane, which is well under control. 
The price which we have to pay is the fact that we have to solve the inhomogeneous BSE 
over the whole range of spacelike instead of the homogeneous BSE at only the on shell 
point = —ml. 

It is the aim of our paper to perform the first analysis of the 2-point correlation function 
in the vector channel (i.e. Tq = 7^) using this approach and to extract the ground state 
mass mg. We will work in the chiral limit (M^ = = Mg = 0), therefore isospin effects 
are not present. Following Ref. we will use a model gluon propagator D{q^) with an IR 
regularized Mandelstam ^ singularity pO| , |3T| at 0. For large an asymptotic free 

UV tail is added. The parameters in D{q^) are fixed in such a way that the weak pion 
decay constant in the chiral limit f^^ = 88MeV is reproduced and, as it was shown in Ref. 
P JTSlJT^ a satisfactorily description of low energy chiral physics and vacuum condensates 
can be achieved. We demonstrate that the dressed quarks which are determined from the 
rainbow DSE are confined in the sense that the quark propagator does not show an 
exponential falloff ~ e~^^^ for large T. We then demonstrate the existence of a bound state 
in the vector channel by showing the exponential falloff of the correlator 11* ^(T) for large T 
from which we can determine the ground state mass mo. 

Our paper is organized as follows: In section 2 we briefly review the definition of the 
dressed quark propagator G and the vector vertex F^ as well as the rainbow DSE for G and 
the inhomogeneous BSE for F^ and demonstrate the validity of the vector Ward-Takahashi 
identity. We then show in section 3 how to obtain the expression for the vector 2-point 
correlator and its Fourier transform in Euclidean space. The numerical results for gluon and 
quark propagators, vertex function and correlator are presented, analyzed and discussed in 
section 4. Finally, conclusions are offered in section 5. 



II. DRESSED QUARK PROPAGATOR AND VERTEX FUNCTION 



A. Rainbow Dyson-Schwinger and Ladder Bethe Salpeter Equation 



Following Ref. ||rT| the inverse of the dressed quark propagator G(p) has the form 

G-\p) = Z7 ■ p[A{p^) - 1] + = Z7 ■ pA{p^) + B{p^) (11) 

in momentum space. The quark self energy dressing S(p) comprises the concept of con- 
stituent quarks by providing a running mass MijP') = B[p^) / A[p'^). It is determined as the 
solution of the rainbow DSE 



S(p) = -g. 



4 2 d^q 



(27r)' 



Dip - qhuGiq)-fu 



(12) 



which is schematically shown in Fig.|I|. For convenience the Feynman type gauge for the 
gluon propagator 



= S.J'^'Diq') 



(13) 
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has been used here, which defines the phenomenological function, D{q^), in Eq.(0). In the 
rainbow DS-BS approximation used in the present work, this function D{q^) for the dressed 
gluon propagator is our main physical input. 

In terms of the components A and B Eg . ([T^) reads: 



A{p)-1 p = Qs - -j—-iD[p - q) (14a) 



In momentum space the dressed vector vertex is determined by the inhomogeneous ladder 
BSE which reads in momentum space 

r,(P, q) = - ^g/ J ^^(^ " K^MK + |)r,(ir, q)G{K - |)7. (15) 

and is schematically displayed in Fig.|[ 

As has been shown in Ref. both the rainbow DSE ( p!2[ ) and the ladder BSE (|I5|) 
can be consistently derived from the action of the global color model (GCM) in an external 
gauge field A/^i^z) using standard functional integration techniques. The main steps of this 
derivation are reviewed in Appendix A and B. 

A crucial consequence is that this formalism ensures by construction invariance under 
local U{1) gauge transformations. At the level of the dressed vertex this invariance is 
refiected in the validity of the Ward-Takahashi identity (WTI): 

g,r,(P, q) = G-\P - |) - G-\P + |) , (16) 



which can be easily verified by substituting it into (|I5|) and using the DSE (|12|) for the quark 
self energy S as well as the definition (|ll|). 

Expanding the r.h.s. of (p!6| ) in q^ and taking the limit q^ ^ ^ leads to the Ward identity 
(WI): 

= ■ (17) 



B. General Form of the Dressed Vertex r^j(P, g) 

Following Ref. |0 in the Feynman-type gauge (0) for the model gluon propagator the 
most general form for the vertex function F^(P, q) which fulfills (|T5|) reads 

F^(P, q) = 11A(,1)(P, q) + «7.A(5(P, q) + t^,^,e^,^f,P^qf,r'^^'\P, q) , (18) 

where for convenience an arbitrary mass scale rj has been introduced on order to render all 
the A*-*-* dimensionless. Furthermore it turns out to be convenient to perform a decomposition 
into longitudinal and transversal components with respect to the external momentum q^: 
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A(^n^,g) = ^^Af + SAf (19a) 



7] 7] 



A(^)(P,g) = Ar. (19c) 



The eight scalar dimensionless coefficients Af (i = 1, . . . , 3) and Af (i = 1, . . . , 5) depend on 
P^, q^ and Cp^, where Cpg = is the direction cosine between P and q. Pj = P^ — 
is the vector transverse to g^, i.e. = 0. The advantage of the decomposition (^) hes 

in the fact that the longitudinal components Af (z = 1, . . . , 3) are determined automatically 
from the quark propagator G by means of the WTI (|T^) 

^^\[ = B(P^) - BiPl) (20a) 

7] 

^A^ = A(P^) - A{Pl) (20b) 
2A^ = A(P!) + A{Pl) , (20c) 

where P± = P± |. 

This leaves only the 5 independent transversal components Af (z = 1, . . . , 5) to be deter- 
mined as solutions of the inhomogeneous BSE (|T^). In order to do this one has to project 
out the single Af for each i which can be done by multiplying (|l^) and (|l^) with appropriate 



Dirac matrices and taking traces. Performing the same operations on the r.h.s of (P^ leads 



finally to a set of 5 coupled linear inhomogeneous integral equations for the A/ (i = 1, . . . , 5) 



X,=U^,X';{P\q\Cl^) + 

POO P + l 

+ dKK' j ^ dCK,M^M^ K\ Cl^)X^{K\ g^ Cl^) . (21) 

The indices i and j are both running from 1 to 5. The vector and the nonzero components 
of the matrices Mij and M.ij have been derived in Ref. |]^. For completeness and because 
we need them later for the determination of the correlator they are also listed in Appendix 
C. 



C. Asymptotic Behavior 

For large values of the quark momentum P^ ^ 00 or for large values of the external 
momentum q^ ^ 00 the propagators reach their asymptotic free forms and therefore the 
kernel in (pTS]) vanishes, which implies that also approaches the naked vertex. This means: 
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or 



and 



or 



together with y4(P^ 





r^(P,g)|p2^ 


oo = 


(-07m 




oo, g^, Cpg) 


= 


, z = l,2 




— s> oo, g^, C|,q) 


= 1 


,2 = 3 




2 ^2 ^ 

OO, g , OpQ^ 


= 


, z = 1,2,4,5 


Af(p2 


2 ^2 ^ 

^ OO, g , o 


= 1 


,2 = 3, 




r^(i^,g)l.^^ 


oo ' 






g^ ^ oo, Cpq) 


= u 


, z = 1,2 




— i 


,1 = 6 




g^ ^ oo, CIq] 


= 


, 2 = 1,2,4,5 




oo, CIq] 


= 1 


, z = 3 


oo) = 


1 and P(P2 - 


> oo) 


= 0. 



(22) 

(23a) 
(23b) 
(23c) 
(23d) 

(24) 



(25a) 
(25b) 
(25c) 
(25d) 



D. Constraints on Vertex Function 



The nonperturbative part of the vector vertex function from Eqs. (0,^ can be constrained 
by the work using QCD sum rules with nonlocal condensates. T^{yi,y2; z)'^^ , defined as 
the second term in Eg. (P8|) has also been shown to be given by the four-quark nonlocal 



condensate 27,28 



Tf,{yi,y2;z)^^ = (0| : q{yi)q{z)j^q{z)q{y2 
Defining the nonlocal quark condensate by 

{0\ : m<liy) ■■ \0) = g{y') (0| : g(0)g(0) : |0) 



|0) 



(26) 



(27) 



the function (7(2/^) gives the space-time structure of the nonlocal condensates. In particular, 
the vertex function r^(0, 0; z)^^ , needed in the three-point determination of the vector 

7!\) with the approximation of vacuum factorization 



vacuum susceptibility is given by Eqs. 

as 



r^(0,0;2;; 
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d'yg{y')9{{z - yf 



(28) 



The vacuum factorization approximation has been found to be satisfactory for deriving 
vacuum susceptibilities [^^. Using the forms of giy"^) such as those studied in Refs. 
1^,|27|,E8|1 one obtains constraints on the dressed vertex functions and thus further constraints 



on the gluon propagator. This will be the subject of future work. 
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III. VECTOR CORRELATOR 



A. Definition and Properties 

We are now ready to study the main object of our investigation which is the correlator 
of 2 vector currents j^{x) = q{x)'y^q{x). In momentum space it is defined as: 

n^,^, (q) ^ J d'xe'^^ {Tj,, {x)j,, (0)) . (29) 

It can be expressed in terms of the dressed quark propagator G and the dressed vertex 
function r^(P, g) as 

n.M(?) = ^ / -^J^ [^,iP^ g)G(P+)7.G(P_)] , (30) 

which is schematically displayed in Fig.^. The trace Tr goes over Dirac and color indices. 
A derivation of (|30|) from the GCM is given in Appendix D. It should be stressed that (|30D 
is a mean field result. 



The WTI (16) for the dressed vertex ensures that the correlator n,y^(g) is transversal 

g^n,^(g) = (31) 

or 

n.,(g) = {s,^ - ^) U{q') . (32) 

This can be easily verified by putting (^) into (^0]) and reflects again manifest invariance 
of our approach under U{1) gauge transformation and therefore conservation of the vector 
current j^{x). Because of ( |52D it is sufficient to consider the contraction n^^(g^). 

B. Spectral Representation and Bound States 

The Kallen-Lehmann spectral representation for 11^^ (g^) is obtained by saturating 
the correlator ( ^9]) with a complete set of on-shell mesonic bound states \I{mj, kj, Xj)) 
(/ = 0,1,2,...) with mass mj, 4-momentum Pj [Pf = —m'j in Minkowski space) and 
polarization Xj: 

poo 'Y 

^M) = l dsp,,{s)j—^^ (33) 

where 

Pi^M) = Pmm(?^)^(?o) 

= (27r)3^5W(P,-g)(0|j^(0)|/(m,,P,,A,))(/(m,,P,,A,)|j^(0)|0) . (34) 

The coupling between the on-shell state |/) and the interpolating quark current j'^ is 
denoted by fj and defined by 
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{Q\j,{x)\I{mi,Pi,\i)) = fie'^^'e^iPiAi) , , (35) 

^(27r)3(2^P,) 



where €^{Pi, A/) denotes the polarization four vector and (cjp/)^ = + kj. Inserting ( pSf ) 
in (|3^) leads in Euclidean space to 

= 3 E • (36) 

In order to filter out the contribution from the ground state |/ = 0) we take the Fourier 
transform of (^) with respect to the Euclidean time T 

2 



and consider the limit T ^ oo 



n^/(T^oo)^^^e-™°^. (38) 



C. Determination of H^^ in Terms of G and 



The final task which remains in order to evaluate 11^^ is to write it in terms of the 
quark self energy functions A and B (0) and vertex functions Xj This is done 

straightforwardly by inserting the definitions (|ll],|18|,|19D into (|30D and carrying out the matrix 
trace. The result is 

n,,(g2) = (-AT,)^ / dPP' / dCp,x ■ { - — xVA^P^ g^ cl^) 



7]^ 



p2 



p2 



-[m-2^—x'T]\l{P\q\C' 



Pq) 



p4 2 



-X 



77° 



CpqTX^ {P ,q , Cpg 



~2^x^TXl{P\q\Cl^)] 



(39) 



where the abbreviations from Appendix C have been used. 

At this point we have to say how to handle the UV divergences in 11^^ (g^). For large 
quark loop momenta P the quark propagator and vertex functions reach their asymptotic 
forms specified in section II C and the loop in Fig.^ is the free loop which is logarithmi- 
cally divergent. In perturbation theory it can be renormalized in the standard way using 
dimensional regularization |3^. In our case it is actually not necessary to evaluate 11^^ (g^), 
because we need only its Fourier transform n^^j*(T), which is UV finite, as any correlator in 
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coordinate space is UV finite. The reason for tfiat is tfiat tlie quark propagator, wliicli falls 
off like p in momentum space, will fall off like e~^^ after Fourier transform, and therefore 
the momentum integral dPP^ in (|39| ) is UV finite. When evaluating 11^^* (T) from ( |39[ ) 
it is therefore convenient to perform the Fourier transform ^^mT ^^fore the momen- 
tum integral dPP^, which avoids a conceptual and technical intricate renormalization 
procedure. 

IV. RESULTS AND DISCUSSION 

Our model gluon propagator has the form 

D{q') = Djniq') + Duviq') , (40) 



where 



gs'Djniq') = (4^'rf)^I^ (41) 

9s'Duv{q') = ^ . (42) 

g2 + LnU#i- + r 



2 



QCD 

The first term Djji (^TD, which dominates for small g^, is a regularized Mandelstam ansatz 
||30| , |3l|| with a strength and an IR regulator A, which models the IR strength of the quark- 
quark interaction. The second term Djjv (|42|), which dominates for large is an asymptotic 
UV tail and matches the known 1-loop renormalization group result with d = ggj'^gAfj ~ M ' 
Aqcd = 200MeV and r = e 0. The model parameters x and A are adjusted to reproduce 
the weak pion decay constant in the chiral limit /^r = 88MeV. We are using 3 different 
parameter sets: 

(43a) 
(43b) 
(43c) 

The forms of D{q^) are displayed in Fig.^, which also shows for comparison the pure UV 
form Duviq"^)- As one can see in all 3 cases the asymptotic form is reached already for 
0.2GeV^ In Ref. § it has been shown that with those values a satisfactorily description 
of all low energy chiral observables can be achieved. 

The second step is the determination of the dressed quark propagator ([TI|) as solution of 
the coupled set of integral equations ( p.4a| ) and ( |14b| ). Fig.^ displays the running constituent 



Set 1: 


A 


= 1.0* 10^ 


-^GeV^ ; 


X = 


1.02GeV 


Set 2: 


A 


= 5.0* 10" 


-^GeV^ ; 


X = 


1.12GeV 


Set 3: 


A 


= 1.0* 10" 


-=^GeV^ ; 


X = 


l.lSGeV 



mass M{p'^) = at spacelike > 0. In order to demonstrate that the quarks are confined 
one has to show that the quark propagator G{p^) has no poles in the timelike region < 0. 
This can again be done by studying the Fourier transform with respect to the Euclidean 
time T 0. It is sufficient to consider the scalar part Gs'- 

aim = r = r '^e'^^' , . (44) 



27r ' J-00 27T q^A^{q^) + B^{q^ 
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If G{p) had a pole at = — the Fourier transform G*{T) would fall off as e for 
large T or 

Ln[G:(r ^ oo)] ~ -MT (45) 

As we can see from Fig.^ (— )Ln[G*(T)] is far from rising linearly with T, at large T, it 
actually rises quadratically with T. This shows that the the parameter set ( ^5D produces 
indeed confined quarks. 

The third step is now to solve the inhomogeneous BSE, i.e. the coupled set of 5 in- 
homogeneous linear integral equations (pT]). We follow the numerical procedure from Ref. 
0. One applies a Gauss quadrature using about tik ~ 30 Gauss points for dK and 
ricjf ~ ^Ckt ~ 10 Gauss points for both J^j^ dCxq and dCxT- It should be noticed 
that the angular integration over Ckt is running only over the the gluon propagator D and 
not over the functions \i{K'^, g^, C]^q). This transforms the integral equations into an high 
dimensional set of ordinary inhomogeneous linear equations for the Aj at the Gauss points, 
which is solved by inverting the (5 ■ hk ■ ncj^^) x (5 ■ uk ■ ncj^^) ~ 1500 x 1500 coefficient 
matrix using standard packages p3| . 

Other than in Ref. 0], where the quark self propagator G which appears in the kernel 
of the BSE (|T5]) has been parameterized by a simple analytic form without solving the DSE 
(0) explicitly, we are using the numerical solutions A and B of (p!4|) in the kernel for the 
BSE ([T5|). Our treatment is therefore a fully consistent solution of both the rainbow DSE for 
the quark propagator G and the inhomogeneous ladder BSE for the vector vertex function 

using the same model gluon propagator D specified in (^) as input. As one can see 
from (|l^) and (|21|) this involves a numerical interpolation of A{P'^) and B^P"^) to the points 
P±. Fig.|^ displays the coefficient functions Aj(P^) for the parameter set 1 ( [43 a| ) for a fixed 
spacelike value of = O.lGeV^ and a fixed angle Gpg = 0. As one can see our result is 
qualitatively similar though quantitatively noticeably different than the one from Ref. [|^. 
Similar than Ref. [0 we also find that the Aj(P^, g^, Gpq) are only very weakly dependent 
on the angle Gpq if P^ and are fixed. 

Finally Fig.^ shows our main result, the Fourier transformed correlator 11*^^ (T) for large 
Euclidean times T evaluated with the three sets (|43|). As mentioned earlier it is essential to 
evaluate the Fourier integral ^^mT ^gfore doing the integral over the quark momentum 
dPP^ in order to avoid UV infinities. From Fig.H we see clearly that 11*^^ (T) has 
an exponential falloff ~ e~'^''^for large T and moreover that the ground state mass mo is 
practically the same for all three parameter sets, because the curves are parallel. The ground 
state mass mo can be extracted by carrying out a linear fit to the numerical result. One 
finds mo ~ 590MeV for all three sets. This has to be compared with the phenomenological 
values of 780MeV for the ground state vector meson mass. The instanton liquid model 
finds mo = 950MeV, whereas quenched lattice QCD |^ obtains mo = 720MeV. 

V. CONCLUSIONS 

To summarize: We have studied a qq bound state in the vector channel using a nonlocal 
and confining model quark-quark interaction, which respects all global symmetries of QCD. 
The model parameters have been chosen to give a good description of low energy chiral 
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physics for the Goldstone degrees of freedom (vr, rf). For the first time we have demon- 
strated the existence of a higher mass state by evaluating the correlator of 2 interpolating 
vector currents in Euclidean space and showing the exponential falloff of its Fourier trans- 
form at large Euclidean times. This employs a consistent treatment of the dressed quark 
propagator G and the dressed vector vertex F^, which are both determined from the model 
quark-quark interaction by the rainbow Dyson-Schwinger equation for G and the inhomoge- 
neous ladder Bethe-Salpeter equation for F^. The method stays until the end in Euclidean 
space and all momenta are spacelike. This avoids any unsafe analytic continuations of the 
dressed propagators into the timelike region. We have found that the ground state mass 
is practically independent on the parameters of the model interaction and lies at about 
590MeV. 

For the creation of a bound state it is obviously necessary to have an highly nonlocal 
quark-quark interaction, which is reflected in its large nonperturbative strength in the JR. 
This leads not only to dynamical chiral symmetry breaking but also to a strong momentum 
dependence of the running dynamical quark mass M(p^) and, in our case even to confined 
quarks. This is an important advantage over with the Nambu-Jona-Lasinio (NJL) inter- 
action which can be regarded as a special case of the GCM with a local quark-quark 
interaction D'^^{x,y) = S'^^5^^5^^\x — y) and has been subject to extensive studies The 
NJL model has a momentum independent constituent mass M and is therefore not confining. 
Though it gives a satisfactorily description of the physics of associated with the Goldstone 
degrees of freedom (vr, K, rf) it has trouble to deal with higher qq bound states, which are 
lying around the 2M threshold and are therefore unstable or barely stable against decay 
into a quark-antiquark pair. 

As long as ground state properties are considered, the method which we have used can 
also be applied to other channels, such as e.g. the axial vector mesons, as well as to 3-point 
functions, which are the basis for the study of meson form factors p6[|. The only implicit 
assumption which we have to make is the existence of a sufficiently large gap between the 
ground state and the first excited state. If ground state and first excited state are too 
close together, the exponential falloff is only evident at very large numerical values of the 
Euclidean time T, where the result is likely to be swamped by numerical errors in the Fourier 
transform due to the highly oscillating behavior of e*"^"*"^. 
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FIGURES 



FIG. 1. Rainbow Dyson- Schwinger equation for the quark self energy S. 



FIG. 2. 2-point correlator 11,^^ in the vector channel. 



FIG. 3. Inhomogeneous ladder Bethe-Salpeter equation for the Vector Vertex Function F^. 

FIG. 4. Model gluon-2 point functions D{s) (]40| ) for three the 3 different parameter sets (^3|). 
In each case the value of the pion decay constant in the chiral limit /^r = 88MeV is reproduced. 
For comparison we also show the form of the UV tail Djjvis) (^2|) (dashed-dotted). 



FIG. 5. The dynamical (constituent) quark masses M{s) for the three parameter sets (43) 



FIG. 6. The Euclidean time Fourier transform Gs*{T) of the scalar part of the quark propagator 
Gg (arbitrary scale) for the three parameter sets (^3|). 



FIG. 7. The five components Ai(-P^), i = 1 ... 5 of the vertex function r^{P,q) as defined in 
([T^ for a fixed spacelike value of the off shell momentum = +0.1GeV^ and angle Cpq = using 
the parameter set 1 ( |43a| ). 

FIG. 8. The Euclidean time Fourier transform n^^*(r) of the vector correlator (arbitrary 
scale) for the three parameter sets (^). Our numerical results are denoted by circles, squares and 
diamonds, respectively. The straight full, dotted and dashed straight lines are the best linear fits. 
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APPENDIX A: THE GLOBAL COLOR MODEL (GCM) IN AN EXTERNAL 

VECTOR FIELD 



Following Ref. |TT| we consider the Euclidean action of the Global Color Model (GCM) 



in an external vector field A^{x) 

Sgcm% q;A] = J d'xd'y |g(x)[7 ■ - ^7.A(x)]g(x) + ^^-(a;)^^ (x - (Al) 
where denotes the color octet vector current: 



= q{x)^J^q{x) . (A2) 



For convenience we will employ the Feynman-type gauge ([T3|) [0,^ for the model gluon prop- 
agator D. Applying the standard bosonization procedure Il36|j37| the generating functional 



Z[A] = e-^[^l = J PgPge-^«^«['^''''^l (A3) 

can be rewritten in terms of the bilocal auxiliary fields B^{x,y) 

Z[A] = I VE'e-'^^ff^^"'^^ (A4) 

with the effective bosonic action 

S.fAB'; A] = (-)TrLn^-MS^ A] + f dVy^^^^^^f^ (A5) 

J 2g/D{x - y) 

and the quark operator 

g-'[B'; A] = [^-d,- z7.A(x)] S{x -y)+ A'B\x, y) . (A6) 

The matrices A^ arise from Fierz reordering the current- current interaction in (ST) and are 
given by 

A» = i (l,. r,,, ^7.7.) <8 (^1.. ^Aj) » (lie. ^A?) . (A7) 

In the mean field approximation, which is the leading order in the fields B^{x,y) are 
substituted by their vacuum values B^^x^y) which are given as the stationary points of the 
effective action 



SS, 



eff 



SB 



= (A8) 

So 



or 

B',[A]ix,y) = gs'D{x - y)tT[A'go[A]-\x,y)] , (A9) 
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where ^{x,y) denotes the inverse propagator with the self energy S(x,?/) = A^BQ{x,y) in 
the external background field A{x) 

go[A]-\x, y) = [^.d,- tiuA,{x)] 6{x -y)+ A^i3^(a;, y) . (AlO) 

We want to stress that both Sq(x,?/) as well as QQ^(x,y) have an implicit dependence on 
the external background field A{x). If the external field A is switched off Qo goes into the 
dressed quark propagator G = Qo[A = 0] which has the decomposition (p]) : 



G-\p) = i-1 ■ v\A{j?) - 1] + = Z7 ■ + B{p^) (All) 

in momentum space. The quark self energy dressing is determined as the solution of 
the rainbow DSE for ^ = 0) rendering (12): 

S(p) = \9s^^.^^V - q)luG{q)^. . (A12) 



APPENDIX B: VERTEX DRESSING AND INHOMOGENEOUS 
BETHE-SALPETER EQUATION 

In coordinate space the dressed vector vertex T^{yi,y2; z) is given as the functional 
derivative of the inverse quark propagator ^of-^]^^ ( |A1Q| ) with respect to the external field 
Afi. 



SGo[AY\yi,y2) 
5A^{z) 



(Bl) 



Taking the functional derivative in ( AlO ) gives for (|E]] 



T^{yi,y2\z) = {-i)-f^S{yi - y2)5{yi - z) + 



SE[A]{yuy2) 
5Aiz) 



(B2) 



The second term on the r.h.s. of ( |B2|) can be determined by employing the stationary 
condition ([A9D, which after Fierz reordering can be cast into 



(2/1,2/2) = -gs^D{yi - y2)-fu 



SGo[A]{yi,y2) 



SA^,{z) 



(B3) 



A=0 



In order to find an expression for 



5A^(z) 



A=0 



in terms of the quark propagator G. we 



write the definition ( |A1U| ) schematically as 

G^^[A] = G-^ + A^V^ 
which leads to the formal expansion 

go[A]=G-GA^V^G±... . 



(B4) 



(B5) 
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Putting (PD into (RSp gives 



5A{z) 



J A=0 



4 [ 

--gs^D{yi-y2) J duidu2-fuG{yi,Ui)T^{ui,U2;z)G{u2,y2hu, (B6) 



which, after substituting the result into ( p2D renders the inhomogeneous BSE for r^(?/i, ^2; 2;) 
in coordinate space 



^i,{yi,y2;z) 



(B7) 



(-^)7m'^(i/i - y2)S{yi -z) - -Qs D{yi - 7/2) J duidu2'yuG{yi, Mi)r^(ui, U2; z)G{u2, y2)lv 
Fourier transform leads then to the momentum space form (p!5|). 



APPENDIX C: MATRIX ELEMENTS OF THE BS KERNEL 



In this appendix we list the explicit forms for the quantities which appear as kernel of 
the 5 transversal inhomogeneous BS integral and will be also needed in (|39|): 



X = (0,3,1,0,0) 



(CI) 



Ml = 1 ^22 = - 
Ms =3 ^32 = - 
7V33 = 1 ^42 = 1 



(pT)2 



AA44 = M 



'55 



1 



(C2) 



and 



-Mil 

M21 
-M23 
-M25 
-M32 
-M34 
-M41 
-M43 
-M53 



2tx D^F^fr 
-2nDi^V 



7rD2^ {f, - 2^T 
-27rD2^CLT 

TcDl pLpT C 

—2itDi pLp'l 'CKqT 
K 



4^2 



We have used the abbreviations: 
F,^{i-K') ^T + U 



M12 
Mu 
M22 

M2i 

-M33 
-M35 

Mi2 
M55 



2i.b,^,(y+^^ci^w\ 

nDo^x' {f, - 2f T 
-27rDo^x2Ci,T 
~7rD2Kv 



-TT (F,Do - 2^TD2 
TT (-Dq^; 
—TcDi pLpT GKqF2 
—TtDi pLpT Gxq^Fs 

nD^^Fo 



^0 
F2 



i^i + 2 



U E 



Fi^ 



hT + U 



2{i-K^Cl^)^ 



a{K+)l3{K-)-<K-)l3{K+) 
2Kq 



(C3) 



(C4) 
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as well as 



(C5) 



and 



l[ d^KTj^ gs'D + K^- 2PK^1 - ClfKT - 2PKCp,Ck,) {CktT ■ 



The Ckt is the direction cosine between and and 



x^Jl- C' , Pl ^ PCpQ , Pt = PJl - C\ 



'Kq ■ 



(C6) 



(C7) 



APPENDIX D: 



In this appendix we give a derivation of Eq.(^). The time ordered product of the two 
vector currents can be formally written as functional derivative of the generating functional 
Z[A] ( [A3| , [A4| ) with respect to the external vector field A: 



5Af,^{zi)6A^.,{z2 



[Dl] 



J^=o 



At the mean field level the integration in over all possible configurations of the bosonic 
auxiliary field B is substituted by the stationary configuration Bq and we can therefore 
write: 



Z[A] 



where 



Wo[A] = (-)TrLn^o[^]"' + / d'^xd^y 



Bl{x,y)Bl{y,x) 
2gsW{x-y) 
This implies for the functional derivatives 

m^] _ ._jm\Al,-mA] 



5Au{z 

svMA 

5AJz) 



+ 



mo[A] 

6Bo[A] 



5Ay{z) 
' >9>Vo[^] 

dA,{z) 
' dWo[A\ 

dAu{z) 

+i){z\ii [go[A]-i,] \z) 



5B^[A] 
5AJz) 



mo[A] 

5Ay{z) 

' 5Z[A] 
5Au{z) 

5A^^{zi)5A^^{z2) 







J^=o 



5(2) Wo [^] 



5A^,^{zi)5A^.Xz2) 



(D2) 
(D3) 



(D4) 



AA=0 
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Therefore we find for the correlator using 

{Ti^,^{zi)i^^{z2)) =iNc j d'^yid'^y2tT^[-i^^G{zi,yi)V^,^{yi,y2]Z2)G{y2,zi)] , (D5) 



which, after Fourier transform, gives (|30[) . 
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